Evidence of electron fractionalization in the Hall coefficient at Mott criticality 
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Hall coefficient implies the mechanism for reconstruction of a Fermi surface, distinguishing com- 
peting scenarios for Mott criticality such as electron fractionalization, dynamical mean-field theory, 
and metal-insulator transition driven by symmetry breaking. We find that electron fractionalization 
leaves a signature for the Hall coefficient at Mott criticality in two dimensions, a unique feature 
differentiated from other theories. We evaluate the Hall coefficient based on the quantum Boltzman 
equation approach, guaranteeing gauge invariance in both longitudinal and transverse transport 
coefficients. 
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I. INTRODUCTION 

Since reconstruction of Fermi surfaces has been dis- 
cussed in the context of high Tc cuprates it is 
considered to play an essential role in understanding the 
mechanism of metal-insulator transitions J4L Its infor- 
mation is reflected in the Hall coefficient 043' expected 
to distinguish competing scenarios for Mott criticality 
based on electron fractionalization [8], dynamical mean- 
field theory (DMFT) [1], and metal-insulator transition 
driven by symmetry breaking (SB-MIT), respectively. 
Indeed, the Fermi-surface reconfiguration has been re- 
garded as one of the central problems in heavy-fermion 
quantum criticality [l3] , claimed to differentiate compet- 
ing theories |ll| - [l4j such as the Kondo breakdown sce- 
nario [l6j . the local quantum criticality ansatz [l7| . 
and the spin-density-wave theory [l8| . 

In this paper we propose evidence of electron fraction- 
alization at Mott criticality from the Hall coefficient. An 
essential feature of the fractionalizaion scenario is ap- 
pearance of novel charge fluctuations (holons) near the 
Mott critical point, giving rise to an interesting metallic 
state in two dimensions. Such charge fluctuations turn 
out to contribute to the Hall coefficient additionally on 
top of the Fermi-liquid value given by fermion excita- 
tions (spinous). As a result, we uncover two-step jumps 
in the Hall coefficient at the two dimensional Mott criti- 
cal point. On the other hand, both DMFT and SB-MIT 
allow only one-step jump from the Fermi-liquid value to 
infinity at the Mott critical point. See Fig. 1. 

A scaling theory for fermions has been proposed near 
Mott criticality based on an ansatz of electron fractional- 
ization, where the scaling expression of the spectral func- 
tion is determined by both the dynamical critical expo- 
nent and the anomalous dimension of fermions [l^. Al- 
though this scaling theory serves a convenient framework 
for quantum criticality in thermodynamics and trans- 
port, it is difficult to find such critical exponents in the 
Fermi-surface problem [l^, [2l|, where they are intro- 
duced phenomenologically. This difficulty does not arise 
in the Hall coefficient, not renormalized by strong cor- 
relations, thus suggesting an undisputable evidence for 
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FIG. 1: A schematic diagram for the Hall coefficient. Elec- 
tron fractionalization gives rise to two-step jumps in the Hall 
coefficient at the Mott critical point in two spatial dimen- 
sions (left), where novel charge fluctuations (holons) appear 
to contribute to the Hall coefficient in addition to fermions 
(spinous). On the other hand, either DMFT or SB-MIT al- 
lows only one-step jump (right). FL means Fermi liquid and 
MC denotes Mott critical point. 



electron fractionalization as the mechanism of the Fermi- 
surface reconstruction. 



II. A SIGNATURE OF ELECTRON 
FRACTIONALIZATION IN THE HALL 
COEFFICIENT 



We start from the Hubbard model 



H = -t^{cl„Cj^ + H.c.) u^ni^n,^. (1) 



Cicr is an electron field, and = cl^Cia is its density, t 
is a hopping parameter, and u is an interaction strength 
at site i. In this study we focus on the case of two spatial 
dimensions. 

A way to describe electron fractionalization is to ex- 
tract out collective charge dynamics explicitly from cor- 
related electrons, referred as the slave-rotor representa- 
tion [22I, [2^ . Such charge fluctuations are identified with 
zero sound modes in the case of short range interactions 
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while plasmon modes in the case of long range interac- 
tions. Actually, one can check that the dispersion of the 
rotor variable is the same as that of such collective charge 
excitations. The Mott transition is described by gapping 
of rotor excitations in the slave-rotor theory. 

The slave-rotor theory is based on decomposition of an 
electron field as follows = e"*^'/^^, where e"'**' is a 
holon field to carry an electron charge, and fia is a spinon 
field to do an electron spin. Resorting to this fractional- 
ization scenario, one can rewrite the Hubbard model in 
terms of holons and spinous, interacting via gauge fluc- 
tuations. Gauge-field excitations are associated with col- 
lective spin fiuctuations, differentiated from those in the 
spin-fermion model (l8l |. Performing the continuum ap- 
proximation, we reach the following expression for the 
Mott transition 

2mf 

(2) 

where the last term with the field strength tensor /^^ — 
d^ttu — di^Qfj, describes the Maxwell dynamics of gauge 
fiuctuations, and the gauge-matter coupling constant is 
denoted as g. Tnf(p) is the band mass of spinous (holons). 
A detailed derivation of this procedure can be found in 
Refs. dllll. 

This effective U(l) gauge theory has been discussed in 



the Eliashberg framework, where self-energy corrections 
of spinons, holons, and gauge fields are incorporated self- 
consistently, but vertex corrections are neglected [13, [1^ ■ 
Integrating over both spinons and holons, one finds an 
effective action for gauge fluctuations 

ai{q,iVt){5i.j - '^^^a.j{-q,-iVi), (3) 

where multi-scale quantum criticality is uncovered from 
two dynamical critical exponents, z = 1 from critical dy- 
namics of holons and z = 3 from particle-hole excitations 
near the Fermi surface of spinons. 7 is the damping co- 
efficient, proportional to the density of states of spinons, 
and X is the diamagnetic susceptibility of spinons. ag is 
the holon conductivity. Resorting to this gauge propa- 
gator, one can evaluate self-energy corrections of spinons 
and holons self-consistently in each regime. It was shown 
that the holon self-energy due to gauge fiuctuations is 
not relevant, implying that the Mott transition falls into 
the 3D-XY universality class [lljIlHl- On the other hand, 
the spinon self-energy gives rise to scaling for thermody- 
namics and transport in each regime. 

In order to evaluate both longitudinal and transverse 
transport coefficients, we resort to the quantum Boltz- 
man equation approach |26i], given by 



-{vi + afc3f?EL(fc,^)} • [B X dk)G<{k,uj) 
c ■' 

+ee^{vi + dum:L{k,uj)}Tf{k,uj) + {^-4-m:L{k,uj)}duTii^^^ 

liuik,^) = i[2Tf{k,u)G<{k,u) - ^<ik,uj)Af{k,u)] (4) 
for spinons and 

+eiE + S) ■ [{vl + dk^Y.l,,{k,Lo)}T,{k,u) + {c. - - SRE^^^.lfc, c.)}9fcrfe(fc, c.)] [9„nfc(c.)] [A(fc, a;)]^ = ll,i{k,uj), 
Icouik,^) = t[2rt{k,iu)G<{k,,o) ~ ^<{k,oj)A,{k,iu)] (5) 



for holons, where G^j-j^^ (/c, cj) and E^j.j^j(fc,w) are lesser 
Green's function and self-energy of spinons (holons), re- 
spectively, and Aj(f,) (fc, oj), T f(^ii){k,LL>), and ^T,ll,f\k,LL>) 
are imaginary parts of retarded Green's function, self- 
energy, and the real part of the retarded self-energy, re- 
spectively. I'fc^^^ and are velocity and dispersion 
of spinons (holons), respectively. rij(b)(a;) is the Fermi- 
Dirac (Bose-Einstein) distribution function. E and B are 



applied electric and magnetic fields while £ and 13 are in- 
ternal fields related with fractionalization. Since spinons 
do not carry an electric charge in our assignment, they 
couple to internal fields only in a gauge invariant way. 
Equations (4) and (5) are derived in appendix A. 

Inelastic scattering with critical fiuctuations gives rise 
to the collision term of the right-hand-side, where each 
lesser self-energy is given by 
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dv 



k X q 



5L»a (<7, 1^) [{nb{i^) + l}Gf {k + q,uj + iy)+ nb{iy)Gf {k + q,uj- i^)] , 



k X q 



nib 



a {q, v) + ^Dx {q, i^)^ [{ribiiy) + l}G<{k + q,uj + + nfc(i/)G< (k + q,uj- u)] . 



(6) 



r 



QDa(g, v) is the spectral function of the gauge propaga- 
tor, given by z = 3 in the paramagnetic insulating state 
and z = 2 at the Mott critical point. QD\{q,i') is the 
spectral function of the A propagator, where excitations 
of A correspond to potential fluctuations for holons in the 
nonlinear a model description [23| • Such fluctuations are 
essential, the Mott transition belonging to the XY uni- 
versality class [13, HI- 

These Boltzman equations are simplified further in 
the Eliashberg approximation, where self-energy correc- 
tions do not depend on momentum. Inserting the lesser 
Green's function 

G'}(b)^k,uj) = m^(t)(w)Aj(b)(fc,w) 

into the quantum Boltzman equations (4) and (5) with 
Eq. (6), we find so called vertex distribution functions 
for spinous and holons 



Af{kF,uj) 
Ab{k,uj) = 



2 Tic{kF,uj) 1 + iQfTlj.{kF,uj) 

e T^^{k,uj) Ab{k,u:) 

2 ri'^(fc, uj)l + i{nb+ uJb)Tl{k, uj) 



respectively, where Q. 



(8) 

are internal and external cyclotron frequencies, and 
S — -\- iSy and E — Ex -\- iEy are complex fields. 

{k,uj)] ^ = 2T f^b){k,LL>) is the relaxation rate and 

b'tr ''\k,uj)]~^ is the scattering rate associated with the 
transport time, where the former quantity is obtained 
from the Eliashberg calculation, and the latter is found 
with vertex corrections, thus a gauge invariant quantity. 
A detailed procedure on how to solve quantum Boltzman 
equations is presented in appendix B. 

Inserting these vertex-distribution functions into ex- 
pressions for currents, we obtain 



Jf = CfS^p^F^^F 



Jb = Cbc" 



l + infrUT) 



l + i{nb + ujb)Tl{T) 



{£ + E), (9) 



where p^p and Vp are the density of states and Fermi ve- 
locity of spinous with a positive numerical constant C/, 



and Cb is a positive number given by a function of temper- 
ature and the ratio u/t. Cb is a universal number at the 
Mott critical point, and it vanishes in the paramagnetic 
insulator. 

Internal magnetic fields are shifted as 



Xb 



Xf+Xb 



nb 



- i^b 



Xf 



Xf + Xb 



iOb. (10) 



determined from the gauge invariance [28|. Another 
condition for the gauge invariance is the back-flow con- 
straint, 



J/ + Jfc = 0, 



(11) 



implying that the total internal current should vanish 
[231 ■ Then, we reach the final expression for the Hall 
resistivity 



PxyiT) 



1 Xf 



[cr-xyiT)] 



e'^CbXf+Xb^^ ' e^Cfp^pV^p^ Xf +Xb 



Xb 



ujf. (12) 



It is important to notice that the transport time does 
not appear, implying that renormalization of the Hall 
conductivity due to quantum corrections is cancelled by 
that of the square of the longitudinal conductivity. A 
derivation is shown in appendix C. 

In the Fermi liquid phase {u < Uc) holon condensa- 
tion gives rise to %(, — )■ oo Then, we find pxy{T) « 
— — — r^, which coincides with that of the Fermi liq- 

uid state. On the other hand, Xf ^ Xb results both 
at the Mott critical point and in the paramagnetic in- 
sulator. Then, we obtain p^yiT) « resulting from 
only holons. In the insulating phase we see C;, — 0. See 
appendix D. As a result, we obtain Pxy(T — > 0) — > oo. 
However, it is not the case at the critical point, where 
the critical dynamics of holons is characterized by their 
own conductivity. Cb is given by a constant, leading the 
Hall coefficient at the Mott critical point to differ from 
the Fermi-liquid value. 

Integrating over spinon excitations, we obtain an ef- 
fective field theory for critical holon dynamics with the 
z — 2 gauge dynamics given by Eq. (3) and the z = 1 A 
dynamics D\{q,iil) — , where Kg is associated 

with the com pres sibility of critical holons. As discussed 
in Refs. (23. |25||. the z = 2 gauge dynamics does not 
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affect the holon dynamics, preserving the 3D-XY univer- 
sahty class to allow the universal conductivity [13, HH 
at the two dimensional Mott critical point. This gives a 
universal number to Cb, resulting in the universal jump 
for the Hall coefficient at the Mott critical point. On the 
other hand, if spinous arc in the diffusive regime, the 
gauge dynamics is given by Da{q,ifl) = af\n]+'^„q ^^^^ 
the dc spinon conductivity ct/. Such gauge fluctuations 
will affect the holon dynamics, and the Mott transition 
does not belong to the XY universality class any more. 
Then, the universal conductivity depends on the damp- 
ing coefficient a / , thus Cb becomes a function of disorder 
and the density of states of spinous. The Hall coefficient 
does not become universal. However, it must be differ- 
ent from the Fermi-liquid value, giving rise to two-step 
jumps in the Hall coefficient at the Mott critical point. 

We claim that the (universal) jump in the Hall conduc- 
tivity will not arise from either the DMFT framework or 
SB-MIT theory. As revealed in the quantum Boltzman 
equation approach, the Hall coefficient does not renor- 
malize due to quantum corrections. Recalling that the 
configuration of a Fermi surface is unchanged at the Mott 
critical point in these frameworks, the Hall coefficient of 
the Mott critical point remains the same as the Fermi- 
liquid value, and jumps to an infinite directly in either 
the paramagnetic insulator or an antifcrromagnetic insu- 
lator (Fig. I). 

One may criticize the assumption for the existence of 
the paramagnetic insulating state, referred as a spin liq- 
uid phase [321] . Although metal- insulator transitions are 
usually involved with symmetry breaking, such an insu- 
lating phase has been proposed in the DMFT framework 
and the gauge theory scenario [8, j54]. In particular, 
such a state has been observed by various numerical sim- 
ulations when geometrically frustrated lattices are con- 
sidered to prohibit magnetic ordering [35| . However, the 
stability of the spin liquid state is not an important issue 
in our study. The point is whether electron fractional- 
ization appears or not at the Mott critical point. If it 
occurs, the two-step jumps should be shown in the Hall 
coefficient, independent of the nature of an insulating 
phase. 



III. CONCLUSION 



In this study we discussed how the spin-charge sepa- 
ration scenario near two dimensional Mott criticality can 
be verified from the Hall coefficient. An essential point 
is that a novel metallic state emerges at the Mott criti- 
cal point, described by critical dynamics of holon excita- 
tions. As a result, such charge fluctuations contribute to 
the Hall coefficient additionally on top of the Fermi liq- 
uid value, given by spinon excitations. Interestingly, the 
jump of the Hall coefficient turns out to be universal in 
the clean limit, given by the 3D-XY universality class of 
critical holon dynamics. We claim that the two-step be- 
havior of the Hall coefficient (Fig. 1) can be regarded as 
a fingerprint of electron fractionalization because either 
the dynamical mean-field theory or symmetry-breaking 
metal-insulator transition allows only one-step jump. 

One cautious person may criticize that the behavior 
of the Hall coefficient depicted in the right hand side of 
Fig. 1 and that of the left hand side do not differ from 
each other in the experimental point of view, if the Mott 
transition occurs really just at a critical point instead 
of a critical region. In this respect the behavior of the 
Hall coefficient at finite temperatures may be essential. 
Since the jump in the Hall coefficient is given by novel 
metallicity of holon dynamics, we are suspecting an in- 
teresting scaling behavior in U /T and H/T for the Hall- 
coefficient jump, where C/ is a tuning parameter, here the 
Hubbard- J7, and H is externally applied magnetic field. 
In particular, the scaling exponent will be given by that 
of the 3D-XY universality class in the clean limit. Unfor- 
tunately, we cannot find any previous researches on the 
scaling behavior of the Hall coefficient in the superfluid 
to Mott insulator transition of the bosonic Hubbard-type 
model. This will be an interesting and important future 
work. 

This work was supported by the National Research 
Foundation of Korea (NRF) grant funded by the Korea 
government (MEST) (No. 201 1-0074542). 
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Appendix A: To derive the quantum Boltzman equation in the presence of an external magnetic field 

In appendix A we derive the quantum Boltzman equation. This derivation extends that of Mahan's textbook (26j 
to the situation with an external magnetic field. We start from equations of motion for the lesser Green's function, 

[idti - Ho{ri,pi)]G<{xi,X2) = J dx3[T,t{xi,X3)G<{x3,X2) - T,< {xi, X3)Gi{x3, X2)], 

[-idt2 - Ho{r2,-p2)]G'^{xi,X2) ^ / dx3[Gt{xi,X3)T.<{x3,X2) - G<{xi,X3)T.t{x3,X2)]. (AI) 



The retarded Green's function remains the same as that in equilibrium [2; 



First, we focus on the left-hand-side of these equations. These two equations can be rewritten as follows 

[iidi^ - dt,) - J?i - H2]G<{xuX2) = 
[i{dt,+dt,)-Hi+H2]G<{xi,x2) = 

iJi = Ho{ri,pi), H2 = Ho{r2, -pa). (A2) 

The main point is to derive a gauge invariant expression in the presence of both electric fields and magnetic fields. 
Introducing the center of mass (T, R) and relative {t, r) coordinates 

ti,2=T±*-, ri,2 = ii±^, (A3) 
we rewrite the covariant derivative of the Hamiltonian as follows 

pi + eEyti + X n = (p + eEyT + x R) + Up + eEJ + ^B x r), 

P2 - eE^t2 - ^Bxr2 = -{p + eE^T +^BxR) + ^{P + eEj + -^B x r), (A4) 

where the applied electric field is decomposed into Ey and Eg and pi^2 = ±P + is used according to the coordinate 
transformation. 

Inserting the transformed covariant derivative into the Hamiltonian part of the equation of motion, we rewrite Eq. 
(A2) as follows 

2[idt - -^{p+eEyT+-^B x Rf - -^{P + eE^t + x rf + eE,R]G<{x,X) = 

[idr - — (p + eE^T + —B x R)(P + eEj + —B x r) + eE,r]G< (x, X) = (A5) 
m 2c 2c 

where df^ + dt^ = dr and dt^ — dt^ = 2dt are used according to Eq. (A3). 
Performing the Fourier transformation, we see 

2[0 + eEsR - ^{q + eE^T + ^B x Rf + ^(^a + eE^dn + ^B x dgf]G<{q, ft; R,T) = 

i[dT + eEsdg + -{q + eE,T + -^B x R){dH + eE^dn + -^B x dg)]G<{q, ft; R, T) = .... (A6) 
m Zc Zc 

This expression is not gauge invariant. In order to make this expression invariant for the gauge transformation, we 
perform 

n + eEsR -^uj, q + eEyT + —B X R^k, 

2c 

dR^dR + eEsd^ + ^Bxdk, dr ^ dr + eEA- (A7) 
As a result, we find the gauge invariant expression in this transformed coordinate, given by 

[iJ-ek + ^ (Or + eEd^ + -Bx dkf]G< {k, u; R,T) = 
om c 

i[dT + VkidR + -Bx dk) + eE{dk + Vkd^)]G<ik,uj; R,T) = ... (A8) 

c 

Following Mahan's textbook, the right-hand-side of Eq. (A2) can be expressed as follows 

[to-ek + ^{Or + eEd^ + -Bx dkf]G<{k,uj;R,T) 
am c 

= ^J dze-"'^ J dye-"iy[My,X + ^)G<{z,X-^)-E<{y,X + ^)Gi{z,X-^) 

+G,{y, X + |)E<(z, X - |) - G< (?y, X + |)EK^, X - |)]g/, 
i[dT + Vk{dR + ^B X dk) + eE{dk + Vkdu,)]G<{k,uj;R,T) 

= I dze-'"^! dye-"'y[J:t{y,X+^)G<{z,X-^)-J:<{y,X + ^)Gt{z,X-^) 

-Gt(y, X + |)I]<(^, X - |) + G<(y, X + |)Sf(z, X - |)]g/, (A9) 



where the subscript "GI" represents gauge invariancc, clearer below. 
One can rewrite the self-energy part in a form of the Moyal product 

1= J J dye-'iyi:{y,X + ^)G{z,X - ^) = E{q,X)ei^^^'^^-^^^^^G{q,X), (AlO) 

where the first derivative in the exponential acts to the front part and the second derivative does to the behind part. 

Near equilibrium, variations with respect to positions are not so large, allowing us to expand the exponent of 
derivatives as follows 

/ = E(g,X)G(g,X) + '-[dgE{q, X)dxG{q, X) - dx^{q, X)dgG{q, X)] + ... 

« S(g,X)G(g,X) + '-[i:{q,X),G{q,X)]. (All) 

This is called the gradient expansion. The Poisson bracket is defined as 

X), G{q, X)] = X)dxG{q, X) - dx^{q, X)d,G{q, X) 

= dn^drG - dr^dnG - {dq^dnG - dnJldgG). (A12) 

As discussed in the derivation of the left-hand-side, the above expression is not gauge invariant in the presence of 
external fields. Performing the following transformation for the right-hand-side 

n + eE,R^io, q + eEyT+^B X R^k, dn ^ dn + eE,d^ + x dk, dr ^ dx + eE^dk, {AU) 

we write down the Moyal product in a gauge invariant way 

/ = x')e^^^^^^~^^^^^'^^^^'^^''~^^^"^^^^'^"^^^^^^i^G(k X) 

« E(A:,X)G(A;,X) + ^-{dk'T.{k,X)dxG{k,X) - dx^{k,X)dkG{k,X)] 

+ '-eE[d^E{k,X)dkG{k,X) - dk^{k,X)d^G{k,X)] + ^-^Be,ijdkMk,X)dkjG{k,X). (A14) 
As a result, we find two gauge invariant equations of motion for the lesser Green's function 

[to-ek + ^{dR + eEd^ + -Bx dkf]G< = 3RE,eiG< + S<3RG,et + 7p>,G<] - '-[^<,G>] 
8m c 4 4 

+ '-eE[d^i:>dkG< - afcS<5„G< - d^i:<dkG< + afeE<a„G>] 

+ '-'-B[dk^i:>dk,G< - dk,T.<dk^G< - dk^Y.<dk,G< + dk,T.<dk^G>], 

i[dT + Vk{dR + -Bx dk) + eE{{l - d^m:ret)dk + (vk + dkm:ret)d^} + -Be,,jdk^m:retdk,]G< 
c c 

-ieE[d^E<dk^Gret - dk^<d^mret] " l^Be,ijdk,^< dk,^Gret = S>G< - S<G> + z[3?S^et, G<] + l[E<,mret], 

(A15) 

where the following relations are utilized 

Gt — Gt = 2^Gret, Gt + Gf = G^ + G^, 

St - Sf = 25RS„.f , St + Ef = S< + S> . (A16) 

Near equilibrium, systems are expected to be in homogeneous {OrG'^ — >■ 0) and steady [drG^ — 0) states. Then, 
the Poisson bracket vanishes identically. These two equations can be further simplified as follows 

[w - efc]G< = 5RE^etG< + E<5iG^ef + ^-eE[d^Y> dkG< - dk^<d^G< - d^^<dkG< + 5feE<a„G>] 

+ '-tB[dk^Y>dk,G< - dk,T.<dk^G< - dk^Y.<dk,G< + dk,^<dk^G>], 

-{Vk + dkm:ret){B X dk)G< + eE[{l - d^^ret)dk + {Vk + dk^EreMG< 

c 

-eE[d^E<dkmret - dk^<d^mret] " ^5e,,,5fe. E<5fc^. 3?G,et = E<A - 2rG<, (A17) 



where the following identities of = — iA and = E< — 2iT are used for the second equation. 
The first equation gives rise to the following solution in the linear response regime 



G<{k,uj)=if{L0)A{k,uj) + O{E'';n>l), 



where 



A{k,u!) 



2T 



(u;-efe-5RS,et)2+r2 



(A18) 



(A19) 



is the spectral function with T = —2Q'Eret- 

Inserting this solution into the second equation, we reach the quantum Boltzman equation 

^{vk + dk?fi^ret){B X dk)G< + eE[{vk + 5fcJRS^et)r + {u-ek- JRE„t)5fcr][a„/(w)][A(A:, w)]^ = -i{E<A - 2rG" 



(A20) 



Appendix B: To derive Aj((,)(A:,a;) from the quantum Boltzman equation 

1. To derive Af{k,u)) 
Inserting the ansatz for the lesser Green's function 

Gf{k + q,uj + u) = iAf{k + q,uj + z^)|n/(a; + ( ^ )^fc+g " ^f'^^ + 9,^ + !/)| (Bl) 



into the self-energy expression, we find 



E<{k,w) = iY^ 



dv 




k x q 



^Da{q, v) 



^nh{v) +nf{u) + u)}Af{k + q,u + u) - {ubi-u) + nf{oj - v)}Af{k + q,(^ -v) nf{uj 
k X q 



ruf 



^Da{q, v) {nb{v) + n/(w + v)} "" i^^!^^^^^ ^li^ + q,oj + v)v{^^ ■ Af{k + q,uj + u) 



l-rifjuj-iy) 



where the following identities for thermal factors, 



(B2) 



{ubiu) + l}nf{u! + 1/) = nf{ijj){nb{i') + nf{uj + v)}, nb{i')nf{u) - v) = -nf{u}){nb{-i') + n/(a; - v)} (B3) 



and 



r / N -.-1/ dnfioj + v)\ r / N / m/ dnf{uj)\l — nf((jj + v) 

{n,(.) + = + + -)}(r^) ' 



1 - n/(w) 



(B4) 



are utilized. 

Inserting both the lesser Green's function and lesser self-energy into the quantum Boltzman equation, we obtain 
the following expression 



i-[vl ■ {B X dk)vi] ■ Af{k,Lu) - eS ■vlTf{k,uj)Af{k,uj) = -2Tf{k,Lu)vl ■ Af{k,uj) 

k Q ^ c 

QDa{q,u)Unb{v) + nf{oj + v)}Af{k + q,uj + v)vl^g ■ Af{k + q,uj + i/) 



-{nb{-v) + nf{u} - iy)}Aj:{k + q,u} - iy)vl_^_^ ■ Af{k + q,u!- zy)|. 



(B5) 



where 



dv 



k X q 



|{^^5(^^) + nf{iJ + i')}Af{k + q,ij + u) - {nbi-v) + nf{LO - v)}Af{k + q,LO - u)^ 
is the relaxation rate. 

We write down the above expression as follows, introducing x and y components explicitly 

-i BAfy{k,uj) - [Af{k,ij)]e£x^ f{k,uj) = -2T f{k,uj)Kfx{k,uj) 

nifc 





k X q 


Jo 


ruf 



'^Da{q, {nb{iy) + nf{uj + i^)}Af{k + q,uj + v)^^Kfx{k + q,^^J + v 



fx 



-{nb{-v) + nf{uj - v)}Af{k + q,uj- v)^^hfx{k + q,ij-u)\, 

i BAfx{k,ijj) - [Af{k,uj)]e£yrf{k,ijj) = -2T f{k,uj)Afy{k,uj) 

rrifc 

k X q 



+ V / — %Da{q,u)\{nbiu) + Hficj + u)}Af{k + q,u + iy)^Afy{k + q,ij + u) 

g Jo m/ I. vj^y 

-{nbi-u) + nf{u- i')}Af{k + q,u- v)^^Kfy{k + q,uj-v)\. 

% 

Introducing 

Ax{k,ij) + iAy{k,uj) = A{k,w), £x + i£y = £, 
the quantum Boltzmann equation becomes 

i-^BAf{k,uj) - e£Tf{k,uj)Af{k,uj) = -2T f{k,u:)Af{k,uj) 

VflfC 

2 



r°° dv 


k X q 


Jo 


ruf 



'^Da{q, v)^{nb{v) + nf{u) + i')}Af{k + q,u> + i') cosOAf{k + q,u> + i') 
-{nb{—v) + nf{u) — p)}Af{k + q,u) — u) cos9Af{k + q,uj — u)^, 



where cos 9 is introduced for the velocity ratio. 

Considering the Fermi surface, we perform the following approximation 



( 



1 + i 



eB 



1 



TUfC 2T f{kF, Lo) 



Af{kF,u)) = y U £ + 



1 



kp X qi 



mf 



^Da{q, v) 



2 ■ 2Tj{kF,u:) ^ 

^{nb{y) + + v)}Af{kF + q.,uJ + u) - {nb{-v) + - v)}Af{kF + q,uj - u)^ cos6Af{kF, 
As a result, we reach the final expression for the vertex distribution function 



Af{kF,oj) 



e Tl^{kF,uj) Af{kF,oj) , 
2 ricikF, u))l + iujlTl^(kF,uj) 



where 



^ 2Tf{kF,u:) = dd dcose / —[vi\os^e/2)]QD,{q,u) 

r/c(fcF,a;) ^tt J J-i Jo 

|{nf,(z/) + nf(oj + v)]Af(kF + q, w + z^) - {nb{-v) + n/(a; - v)]Aj(kF + q, w - v)^ 
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is the relaxation rate and 



1 



i^{kF,io) J J-1 JQ 



dv , 



+ nf{u) + v)}Af{kF + q,ui + ly) - {nb{-v) + nf{ui - u)}Af{kF +q,uj - - cos 61] (B13) 



is the scattering rate associated with transport. Note that there is 1 — cos^ factor in the transport time. 
The spinon current is given by 



d?k k,, f du 



resulting in 



(27r)'' mf J 27r 
d'^k f dui / dnf{u}) 



Gf{k,uj), 



UJ), 



r d^k 


f duj j 


dnf(uj)\ 


1 (27r)3j 




duj ) 



Then, we obtain 



fo f d^k f dw . 9^ / dnf(ijj) 
Jf = Jf. + ^Jfy =4'jj^j^ ("^ 



Af{k,u))Af{k,u)) 



^ J (277)3 y 2n "^K duJ ^' TUkF,Uj)l + lL0lTl{kF,Uj) 

rlriT) 



where 



dcosOkp sin^ 9kf 



d^—f 



(B14) 



(B15) 



(B16) 



(B17) 



with a positive numerical constant C. Note that the relaxation time in the denominator is cancelled after the ^ 
integration, making C/ a positive numerical constant, too. The final expression in Eq. (B16) with Eq. (B17) coincides 
with the well known result, justifying our scheme of approximation. 



2. To derive Ab{k,u) 

One can derive the vertex distribution function for holons in a similar way of that for spinous. Inserting the lesser 
Green's function 



G^(fc + g,a; + I/) = iAb{k + q,u) + i/)|n6(w + i/) - ^ 

into the lesser self-energy, we obtain 

r°° / k y. n \ 
j:<{k,u)=iY, ( ^Da{q,u)+QDx{q,iy)) 

.1 (\ ^ TTlh ' 



9a; 



Af,(/c + g,w + i/)} (B18) 



- + v)}Ab{k + g, w + i/) - - - + g, w - i/) n(,(w) 



\nb\y) - nb[uj + '^)|— j— i-^Ab[k + q,uj + • Ab[k + q,uj + u) 

L i -|- Tit) \J-^ ) 

-{nbi-v) - nb{uj - v)} + 9, w - z^X+g • Ab{k + q,uj- i/)J 



(B19) 
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where the following identities, 

{mil') + l}nb{uj + !/)= nb{u)){nb{i') - nb{uj + u)}, nb{v)nb{uj - v) = -nb{co){nb{-iy) - mico - u)} (B20) 

and 

dnb{u) + u)\ / dnb{ui)\l + nb{u) + v) 



n,(.)(-^!^^^) =-{n,(-.)-n,(a. -.)}(- ^n.^-h ^ + - .) 



duo 



1 + nb{Lu) 



(B21) 



are utilized. 

Inserting this lesser self-energy and the lesser Green's function into the quantum Boltzman equation, we obtain 



i-vi ■ [{3 + B)x dk{vl ■ Ab{k,oj)}] - e{£ + E) ■ vlTb{k,uj)Ab{k,uj) = -2r6(fc, ' M{k,w) 



+ 



r-( 


k X q 


Jo V 


nib 



{q,iy) + QDx (q,!^)^ ^^{nb{u) - nb{uj + i')}Ab{k + q,u + v)vl_^_g ■ Ab{k + q, uj + v) 
{nb{-u) - nb{uj - u)}Ab{k + q,aj- z^)t^fc+, • Ab{k + g, w - (B22) 



where 



r-( 


k X q 


Jo V 


rub 



a(q,i^) +'^Dx{q,iy)^ 



|{n6(i/) - nb{uj + i')}Ab{k + q,u) + i')- {nb{-v) - nb{uj - u)}Ab{k + q,LJ - 
is the relaxation rate. 

We write down the above expression as follows, introducing x and y components explicitly 

-i — {B + B)kl{k, to) - e{£^ + E^)Tb{k, io)Ab{k, w) = -2Vb{k,i.)Kl{k, w) 
rribc 



(B23) 



+ 



r-( 


k X q 


Jo 1" V 


rUb 



'?sDa{q, u) + QDx{q, i^))[{nbiu) - nb{uj + u)}Ab{k + q,LJ + u)^A^{k + q,uj + u) 



-{nb{-v) - nb{w - iy)}Ab{k + q,uj- z^)-^A^(A; + q,w-i')\, 

i—{B + B)A.l{k, uj) - e{£y + Ey)Tb{k, w)Ab{k, w) = -2Vb{k, uj)Kl{k, w) 

nibC 



+ 



r-( 


k X q 


Jo 71" V 


rub 



^Da{q, v) + QDx{q, '^))\{nb{iy) - nb{co + iy)}Ab{k + q,u + v)^hl{k + q,uj + v) 

by 

-{nb{-u) - nb{uj - u)}Ab{k + q,cj- z/)^A^(fc + g,a; - i/)}. (B24) 



Introducing 

A^{k,ij) + iAl{k,uj) =Ab{k,iuj), E^ + iEy = E, 
the quantum Boltzman equation becomes 



(B25) 



i {B + B)Ab(k,u)) - e{£ + E)rbik,uj)Ab{k,u)) = -2rb{k,u))Ab{k,uj) 

mbC 

QDa{q,u) + QDx{q,v)]l{nb{u) - nb{uj + v)}Ab{k + q,uj + v) cos9Ab{k + q,u) + u) 



-{nb{—i') — nb{u! — v)}Ab{k + q,ui — cos6Ab{k + q,u} — v)^, 



(B26) 



11 



where cos 9 is introduced for the velocity ratio. 

This expression is further simphfied, neglecting q and u in the vertex distribution function, 

i—(B + B)Ab(k,uj) - e(£ + E)rb(k,u))Ab(k,u)) = -2r b(k,u))Ab(k,uj) 
iribC 

+ V( ^'^Da{q,u) +'^Dx{q,iy)'^{{nb{u) - UbiuJ + u)}Ab{k + q,co + ly) 

-{nb{—v) - nb{uj - v)}Ab{k + q,u) — !/)| cos6Ab{k,u)). (B27) 

Then, we reach the final expression 

A //. . .^ _ '^ ^tr(fc.t^) Ab{k,Uj) I rR9S^ 

2T°^{k,u)) l + i{Ctb+(^b)Tt^{k,io) 



where 



|{n6(z/) - nbicj + i')}Ab{k + q,LJ + v) - {nbi-v) - nb{u) - u)}Ab{k + q,u) - u)^ (B29) 



is the relaxation rate and 



f f f dv / k^ 

/ dqq^ dcose —(—cos^{0/2)QDa{q,,y)+QDx{q,iy) 
Tt^r(k,uj) (2Try Jo J-1 Jo TT Vm^ 

|{n6(z^) - nb{LO + u)}Ab{k + q,uj + 1/) - {nbi-u) - nb{uJ - v)}Ab{k + q,(jj - z^)|[l - cos 6*] (B30) 
is the scattering rate associated with transport. 



Appendix C: To derive the longitudinal and Hall conductivities from the current formulae 

The back-flow constraint results in 

e = -, — — E. (ci) 

Then, the electric current is given by 



Jel = E. (C2 

^^^^^^ l+i^.,riAT)+^'^l+i^.,rUT) 

Cooking the above expression as follows 

CfCbpivi\{{T)TUT) 



Jel = e 



{Cfp^prUT) + Cbr^T)) + (Cfpiv{.\{{T)^^uJbrUT) + CbrUT)^WfrUT)) 
{Cfpivi\i{T)+Cbrl{T) - i(C/p^«/V/,(T)-^.;,ri;(T) +C6ri',(T)-^c^/T/,(T))}(£;. +i£;,)(C3) 

^ Xf + Xb Xf + Xb ^ > 
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we obtain both longitudinal and transverse conductivities, 

C/Ci,p^4V,{,(r)r,UT) (c^p^4V/,(T) + C(,r,';(r)) 

^ (rj.._^2 CfC,p^pvi\l{T)Tl{T) 

^xy\-^ ) — 2 2 

{Cfp^prlAT) + C^tUT)) + (CfpWFTUT)^^uj,TUT) + C,rUT)^UfTl{T)) 
(CfpWF\iiT)-^o;trUT)+C,Tl{T)-^u:fTUT)) (C4) 

^ Xf +Xb Xf +Xb ' 

respectively. These conductivities give rise to the Hall resistivity of Eq. (12). 



Appendix D: To evaluate the holon current 

The holon current is given by 

d^k f duj k"^ 



f d'^k f duj k^ . In ( dnb{w)\ 
Inserting the vertex distribution function into the above, we see 

where 

is the holon spectral function. 

In the Eliashberg theory the relaxation time is assumed to be momentum-independent, Tj'g(fc, ui) ~ tJ'c(w), where the 
singular contribution results from the frequency dependence. Accordingly, the transport time is assumed to depend 
on frequency only. 

Resorting to this assumption, we obtain the current formula 



mU 2tt Jo V duj /V(w-a)^ + [Ti'c(w)]-V T^^{uj)l + i{i 

^ e2 f dij ^ dnb{uj)\ f°° , / r M 1 1 \^ 

mf, J 2n\ duj H (u^+^)r^^^(^) V / Vy^ + l/ 1 

'ml J 2n\ du J J, ' [{oJ + ^i)^Tb,{co)}^z^ + l]^l + i{nl 
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where €$ is a positive numerical constant associated with the angular integration, and Cb is given by 



{T+i,)-'{rl{T)} 



+ 



+ 



dz 



(Z + 1)3/2 1 



r;.(T) 



(Z + 1)3/2 



dz- 



(T + M)3/2{ri',(T)}3 



(r + ,z)3/2{r,^^(r)}3 

e(l-{T + f,)-'{TliT)}-') 



}Q[{T + ^^)-'{Tl{T)}-'-l) 



mt 



2 

+ - 



{^^[^.'0(1^)]-'/' - 5(r + M)'/Mc(7^)}e((T + M)-Hr.'c(T^)}-' - l) 

e(i-(T + M)-HTi'c(T)}-^) 



3(T + M)3/2{r,^,(r)}3 



(D5) 



This expression tells us that Cb vanishes as T goes to zero in the paramagnetic Mott insulator because the presence 
of an excitation gap (/i 7^ 0) implies irrelevance of self-energy corrections in the T — > limit. The only contribution 
resulting from the second term dies out in this limit. On the other hand, the first term contributes to Cb at the 
Mott critical point {fx = 0). The self-energy correction due to A fluctuations is given by J d^qq~^\k + q\~^ 

J_^dcos9 dqq'^q~^{k^ + q'^ + 2kqcos6)~^ = /^^ dcos^ dxx{l + + 2xcos0)~-^, causing [Tgi.{T)]~^ = const. 
in one time and two space dimensions. As a result, we obtain a finite value for Cb- However, it turns out to vanish in 
three spatial dimensions, allowing only one-step jump for the Hall coefficient. 
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